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Abstract 



Azimuthal asymmetries in semi-inclusive production of positive (/1+) and negative hadrons (h^) have 
been measured by scattering 160 GeV muons off longitudinally polarised deuterons at CERN. The 
^ ' asymmetries were decomposed in several terms according to their expected modulation in the az- 

^ [ imuthal angle ^ of the outgoing hadron. Each term receives contributions from one or several spin 

O ■ and transverse-momentum-dependent parton distribution and fragmentation functions. The ampli- 

D I tudes of all -modulation terms of the hadron asymmetries integrated over the kinematic variables 

j^^ • are found to be consistent with zero within statistical errors, while the constant terms are nonzero 

and equal for /i+ and h^ within the statistical errors. The dependencies of the ^-modulated terms 
CsJ ' versus the Bjorken momentum fraction x, the hadron fractional momentum z, and the hadron trans- 

P* I verse momentum pj^ were studied. The x dependence of the constant terms for both positive and 

negative hadrons is in agreement with the longitudinal double-spin hadron asymmetries, measured in 
semi-inclusive deep-inelastic scattering. The x dependence of the sin ^-modulation term is less pro- 
nounced than that in the corresponding HERMES data. All other dependencies of the ^-modulation 
amplitudes are consistent with zero within the statistical errors. 
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Azimuthal asymmetries 



1 Introduction 

Starting from the first polarised lepton scattering experiments at SLAC f\] and at CERN by the EMC 121, 
the longitudinal spin structure of the nucleon has been investigated over the past 20 years by the SMC |[3l . 
HERMES 14], CLAS 10 and COMPASS HQ Collaborations. The cross-section asymmetries Ai and 
A J were measured respectively in inclusive Deep Inelastic Scattering (DIS) 

I+N^t+X (1) 

and Semi-Inclusive Deep-Inelastic Scattering (SIDIS) 

I+N^i' + h+X (2) 

of longitudinally polarised leptons {£) off longitudinally polarised nucleons (N). SIDIS, where in addition 
to the scattered lepton a hadron h is detected, gives access to the individual quai^k spin distributions. From 
the measured spin-dependent asymmetries the contributions of quark spins to the spin of nucleons as well 
as the spin quark distribution functions for valence and sea quarks have been determined. 

The hadron transverse momentum leads to a dependence of the SIDIS cross-section on the hadron az- 
imuthal angle <p (Fig. [T^) and to asymmetries in the hadron production. The asymmetries have been 
predicted lISUTOl and elaborated on in a number of theoretical papers (see Refs. |[TTI - [T6l and references 
therein). The azimuthal asymmetries are related to transverse-momentum-dependent Parton Distribu- 
tion Functions (PDF) and polarised and nonpolarised Parton Fragmentation Functions (PFF). They can 
depend on the transverse or longitudinal component of the quark spin. These asymmetries were first 
observed by SMC HH, HERMES lUSl and CLAS 1191. Further studies were performed both with trans- 
versely polarised targets by HERMES (proton) |[20| and COMPASS (deuteron, proton) ||2ll|22l and with 
longitudinally polarised targets by HERMES (proton, deuteron) ||23l|24| and CLAS (proton) |[25l. Some 
of the measured azimuthal asymmetries (e.g. the so-called "Collins asymmetry") for the proton are rather 
large, reaching up to 10% EOl . others however (e.g. the so-called "pretzelosity"), do not exceed a couple 
of percent ll26l and their investigation requires very high statistics. The asymmetries for the deuteron are 
found to be much smaller than those for the proton indicating opposite signs of contributions from u and 
d quarks iMfBl . 

A search for azimuthal asymmetries for unidentified hadrons with longitudinally polarised deuterons 
is presented below. These new data will test in a wide {x,Q^) range the existence of still unobserved 
azimuthal asymmetries connected to several transverse-momentum-dependent PDFs. 

The Paper is organised as follows. In Section 2 a short theoretical overview with the basic formulae is 
given. The analysis method and the data selection are described in Sections 3 and 4, respectively. The 
results are presented in Section 5 and their stability and systematic uncertainty are discussed in Section 6, 
followed by the conclusions in Section 7. 

2 Theoretical framework 

The kinematics of the SIDIS reaction is shown in Fig. [T^. The 4-momenta of the incident and scattered 
muon are, denoted by i and £', respectively. The momentum transfer is given by q = i — i' with Q^ = —q^ 
and dj is the angle of the momentum q of the virtual photon with respect to the beam. The vectors p'' 
and P|| are the hadron momentum and the longitudinal target polarisation, respectively. Their transverse 
components pj and Pj are defined with respect to the virtual photon momentum. The longitudinal 
component \Pi\ = P\\ cos Qy is approximately equal to P\\ due to the smallness of the angle Qy. The small 
transverse component is equal to \Pt\ = P\\ sin(0y) where sin(0y) « 2{Mx/Q)yJ\ —y, M is the nucleon 
mass and y = q- p/p • i is the fraction of the muon energy lost in the laboratory reference frame. The 
angle ^ is the azimuthal angle between the scattering plane and the hadron production plane and 05 is 
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Production plane 








....^ hi(x) 



(a) (b) 

Fig. 1: (a) The kinematics of the SIDIS process shown for antiparallel target polarisation Ph w.rt. the beam mo- 
mentum {(j)s = 71). (b) Squared modulus of the matrix element of the SlDlS reaction l + N^H' + h+X summed 
over the final states X. 

the angle of the target polarisation vector with respect to the lepton scattering plane. The invariant mass 
squared of the virtual photon-nucleon system W^, the Bjorken variable x and the hadron momentum 
fraction z, characterising SIDIS together with Q^ and y, we. defined as W'^ = {p + qY, x = Q^ /2p ■ q, 
z = p- p^ jp-q where p is the 4-momentum of the incident nucleon. 

The general expression of the total differential cross-section for the SIDIS reaction is a linear function of 
the muon beam polarisation Py^ and of the target polarisation components Pi^ and Pj 

da = daoo + P/idaLo + Pl (daoL + P^daLi) + \Pt \ (daor + P^daLr) (3) 

where the first (second) subscript of the partial cross-sections refers to the beam (target) polarisation. 
The asymmetry a(^) for hadron production from a longitudinally polarised target is defined by 

da^^ - da^^ _ daoL + P^^^Oll - tan Qy (daor + Pii^Olt) 



«(0) 



|PL|(da^^+da^^) 



daoo + Puda, 



(4) 



>UOL0 



where =^ or <^ denotes the target polarisation along or opposite to the muon beam direction and -^ 
denotes the beam polarisation, which is always opposite to the beam direction. The partial cross-sections 
daoo and dai,o do not contribute to the numerator of the asymmetry (Eq. ([4])) while daor and da^r ai^e 
suppressed by the small value of iPrl/lPil = tan Qy k, 2{Mx/Q)yJl —y. 

In the parton model (one-photon approximation, handbag-type diagram) the squared modulus of the 
matrix element of the SIDIS reaction is represented by a diagram of the type shown in Fig. [TJ5. The 
chiral-odd transversity' PDF h\ {x) coupled to the chiral-odd Collins PFF H^{z) is given as an example. 
Each of the partial cross-sections in Eq. ^ is characterised by several terms including a convolution of 
PDF and PFF multiplied by a function of the azimuthal angle of the outgoing hadron. Ignoring the pure 
twist-3 ("tilde") fragmentation functions and retaining only terms up to order {M/Q), the contributions 
to Eq. dU) for unpolarised or spin- zero hadron production have the forms 

daoo °^ xfi{x)0Di{z) + exhf{x)0H^{z)cos{2<p) 
-M 



Q 



+^/2eil + e)—x^(h{x)(S>H^{z)+f^{x)'S)Dl{z))cos^, 



daio 
daoL 



-M 



e(x) ® H^ (z) + g^ (x) (g) Di (z) ) sin (j) 



^2e{l-e)-x 

exh^L (x) ® H^ (z) sin {2(j)) 

+ ^J2e(\^z)—x- (hiix) (g) Hi (z) + ft (x) (g) Di (z)) sin , 



(5) 



' In this Paper we follow the Amsterdam notations I12II16I for PDF and PFF. 
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dOLL '^ Vl - s^xgiL (x) ® D 1 (z) 

+V2£(l-e)— x2 (^gj^ix) ®Di(z) +et(x) <E)H^iz)) cos(/), 

daor °<: e{xhi{x)(g)H^{z)sm{(^ + (^s)+xhiT{x)(g)H^{z)sm{3(l>-(j)s)} 

+x/iV (x) ® Di (z) sin(0 - <ps) , 
dOLT °^ A/l-e^xgi7-(x)(g)Di(z)cos(0-05'), 

where the target spin angle is 05 = or n. The ratio of the virtual-photon flux with longitudinal to 
that with transverse polarisation is given by £ ?» 2(1 —y)/{2 — 2y + y^). The symbol (^ represents a 
convolution of a PDF and a PFF weighted by a function of transverse momenta as defined in Eqs. (4.1- 
4.19) of Ref. |[T6l where also the exact expressions of the cross-sections are given. 

Some of the asymmetry modulations arising from Eqs. ^ and © have already been observed in ex- 
periments with either transversely or longitudinally polarised targets. These asymmetries involve the 
sin0 modulation in doot due to the twist-3 PDFs hi and /^ and that in dOoT due to transversity h\ and 
the Sivers PDF f^. The sin(20) modulation in dOot arising from the "worm-gear" PDF hf^ has been 
seen also. Other asymmetries have not yet been observed experimentally. These include the sin(30) 
modulation in daor related to the "pretzelosity" PDF and the cos (p modulations both in da^r due the 
"worm-gear" PDF gir and in dOit due to the twist-3 gj^ and e^. 

The aim of this study is to evaluate the azimuthal asymmetries in hadron production from the longitu- 
dinally polarised target as a manifestation of the quark-spin and transverse-momentum dependent PDFs 
and PFF mentioned above and to investigate the x, z and p^ dependence of the corresponding modulation 
amplitudes. 

3 Analysis 

The experiment was performed in the muon beam M2 at CERN with positive 160 GeV muons. The beam 
is naturally polarised opposite to the muon momentum with an average polarisation P^ = —80%. 



Briefly, the COMPASS setup 11271 is a two-stage forward spectrometer with the world's largest polarised 
target, tracking detectors and particle identification detectors behind each of the two spectrometer mag- 
nets. Various tracking detectors (MICROMEGAS, GEM, Straw, MWPC, DC) provide a precise deter- 
mination of the particle coordinates, while electromagnetic and hadron calorimeters, muon detectors, 
and a RICH provide identification of secondary electrons, muons and hadrons. The fast trigger and data 
acquisition systems provide for high statistics measurements. 

The method of analysis takes advantage of the polarised target arrangement optimised for asymmetry 
measurements. In 2002-2004 the target consisted of an upstream cell ("U") and a downstream cell 
("D") placed along the axis of a 2.5 T solenoidal magnet centred along the beam direction. The target 
material ^LiD was kept at a low temperature in a ^He-^He dilution refrigerator. The material in the cells 
was polarised longitudinally with opposite orientations by Dynamic Nuclear Polarisation. The beam 
traverses both cells and data are taken simultaneously for both polarisations. To minimise remaining 
systematic effects caused by possible time-dependent variations of the acceptance, the polarisation of 
the cells is reversed three times per day by adiabatically inverting the solenoid field. To avoid possible 
systematic acceptance effects resulting from the different solenoid field orientations, after a few weeks 
the same polarisation configuration is realised with inverted magnetic field. 

For the azimuthal-asymmetry studies double ratios Rf of event numbers are used 

Rf{^) = [</#)/<^.(0)] X [N^j{^)/N^ji^)] (6) 

where M , (0) is the number of events in a given bin originating from the target cell t (t = U, D) with 
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the polarisation orientation p {p = +, — ) and the solenoid field orientation f if = +, — ) w.r.t. to the 
beam direction. 

Using Eqs. ([31 [S]) the number of events can be expressed as 

N'pj = C'fWL'pj f (Bo + fii cos + B2 cos 20 + B3 sin (/) + ... ) 



±Ppj{AQ +Ai sin0 + A2 sin(20) + . . . ; 



(7) 



where C'r{<}>) is the acceptance factor, L' .- is the luminosity, and P* , is the absolute value of the averaged 
product of the measured positive or negative target polarisation and the dilution factor^ calculated for the 
cell t 13. The coefficients Sq, Si, ... and Aq, A\, ... characterise the target-spin-independent and the 
target-spin-dependent parts of partial cross-sections contributing to the denominator and numerator in 
Eq. (|4|, respectively. 

Substituting Eq. (|7]) into Eq. ^, one can see that the multiplicative acceptance factors cancel out as well 
as the luminosity factors if the beam muons cross both cells. The ratios Rf{<p) thus depend only on the 
asymmetry a{<p) (Eq. ([4])), which can be expressed (to first order) by 

afi^) = [Rf{^)-l]/iP^j.+P^^f + P"j + P^j). (8) 

Since the asymmetry should not depend on the orientation of the solenoid field, one can expect a^ = «_. 
From the data a small nonzero difference between a+ and a_ was found. This difference may be due 
to nonfactorisable solenoid-field-dependent contributions in Eq. (|7]). However, these contributions have 
different signs, as it was checked by Monte Carlo simulations, and cancel out in the weighted sum 
a{^) = a+(<^) (Ba^{<p). Therefore, this sum - calculated separately for each year of data taking and 
averaged at the end - is used for the final results. 

4 Data selection 

The data selection starts from the full data set of SIDIS events taken in 2002-2004 with Q^ > 1 (GeV/c)^ 
and y > 0.1. For each event a reconstructed vertex with incoming and outgoing muons and one or more 
additional outgoing tracks is required. Applying cuts on the quality of reconstructed tracks and vertices, 
the target volume, the momentum of the incoming muon (140 < |/| < 180 GeV/c), the fraction of the 
muon energy loss (y < 0.9) and the invariant mass of the virtual photon-nucleon system (5 < W < 
18 GeV/c^), about 96 x 10^ of the SIDIS events remain for further analysis. 

The tracks originating from SIDIS events have been identified as hadrons using the information from 
the hadron calorimeters HCALl and HCAL2 |[27l . A track is considered a hadron track if it hits one of 
the calorimeters, the calorimeter has an associated cluster with energy greater then 5 GeV in HCALl or 
greater then 7 GeV in HCAL2, the coordinates of the cluster are compatible with those of the track, and 
the energy of the cluster is compatible with the momentum of the track. The total number of the hadrons 
is about 53 X 10^. 

The asymmetries are evaluated in the restricted kinematic region of x = 0.004-0.7, z = 0.2-0.9, pj = 0. 1- 
1 GeV/c. The lower x-cut value corresponds to a cut on 2^ > 1 (GeV/c)^, while the highest x value is 
determined by the acceptance. The restriction of the energy fraction of hadrons to z > 0.2 corresponds 
to the cm. Feynman variable Xf ^ z— [pj + ni^/ijW'^) > and assures that the hadron comes from 
the cuiTcnt fragmentation region. This cut removes almost half of the hadron tracks. The high-z cut is 



^The dilution factor is given by the ratio of the absorption cross-sections on the deuteron to that of all nuclei entering the 
target cells. It includes a correction for the relative polarisation of deuterons bound in ^Li with respect to free deuterons. It 
also includes the dilution due to radiative events on the deuteron, which is taken into account by the ratio of the one-photon 
exchange cross-section to the total cross-section. 
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Fig. 2: Distribution of tlie SIDIS events vs. Q^ (left) and vs. y (riglit). 
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Fig. 3: Distribution of hadrons from the restricted kinematic region vs. z (left) and vs. pj (right). 

applied to limit the influence of exclusive channels. The low-p'f cut guarantees a good determination of 
the hadron angles while the high-p^ cut corresponds kinematically to the Xf and W cuts. Distributions 
of the SIDIS events, from which the hadrons are selected for asymmetry evaluations, are shown in Fig.|2] 
vs. Q^ and y. The average value of Q^ is 3.26 (GeV/c)^. The distributions of the hadrons vs. z and p'^ 
are shown in Fig. |3] 



5 Results 

The weighted sums of the azimuthal asymmetry a{^) 
have been fitted by the function 



:fl'+(0)©fl:_(^) for negative and positive hadrons 



a{(j)) 



^^const^^sin: 



'sin(/)+a™20^jj^(2(/)) + a™3^sin(30)+a™^^cos0. 



(9) 



The fit parameters a^^*^^ are connected to particular PDFs and PFFs in Eq. ([5]) and can depend on x, z 
and pj. First, the asymmetries a(0) integrated over these kinematic variables have been calculated and 
fitted in 10 bins. In a second step the dependence of a^^*^^ on each of the variables was studied while 
integrating over the other two variables. For this purpose fits of a{<^) have been performed in 10 bins of 
^ for each of the 6 bins in x, z, and pj. In Eq. ^ we have disregarded the contribution of the cos 2^ 
term which could appear from dcToo in the denominator of Eq. dJ]). This amplitude is expected |[28l to not 
exceed 10% and would enter into Eq. ^ with the factor a^°^^^ which is of order of 10^^ (see Table [B. 
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Table 1: Best values of the a(0) fit parameters for positive and negative hadrons from the five- and one-parameter 
fits. 



Fit parameters 
xlO^ 



const 


23 ±17 


40 ±15 


^sixKJ) 


15 ±23 


-30 ±21 


^sm2(j) 


30 ±23 


-24 ±21 


^sin30 


40 ±24 


-10±21 


^COS(|) 


-4 ±24 


38 ±22 


2/n.d.f. 


6.1/5 


1.0/5 



23 ±16 35 ±15 



10.4/9 7.0/9 



This is beyond our accuracy. The same comments are valid for the contribution of the cos (j) and sin 
amplitudes coming from the denominator of Eq. (01). However, their contributions are automatically 
taken into account by the corresponding fit parameters. 

For the integrated sample the parameters characterising the -modulation amplitudes are compatible 
with zero within 1-1.5 standard deviations. The 0-independent parameters a^°'^^^ differ from zero and 
are almost equal for /i^ and /i+ within the statistical eiTors. The fit parameters are given in Table [T] 
The correlation coefficients of the fit parameters do not exceed 15%. Also given are the results for fits 
considering only a constant term, a{(j)) = a^°^^\ The ^-independent parts of the asymmetries appear due 
to the first term of dOiL, which is proportional to the helicity PDF gn = gi convoluted with the PFF of 
unpolarised quarks in an unpolaiised hadron (see Eq. Q). For the isoscalar deuteron this contribution is 
expected to be only weakly dependent on the hadron charge, as is confirmed by the results. Moreover, 
these constants are related to the hadron SIDIS asymmetry A'I ^{x) for deuterons (see below). 

The dependence of the amplitudes of the modulation on the kinematic variables is shown in Figs. |4]-[8l 
The bin sizes are optimised to contain more than 10^ events in each bin. Some points are slightly shifted 
horizontally for clarity. Only statistical errors are shown. Systematic uncertainties are estimated to be 
much smaller then the statistical ones (see Section 6). 

The asymmetry parameter a^™^'^{x) divided by the average muon polarisation and the virtual-photon 
depolarisation factor Dq ?» |P^|Vl — £^ in the corresponding x bin is shown in Fig. |4l By definition, 
the value a^"^^\x)/Do is equal to the asymmetry A^ ^^(x) published earlier |f7l|. The agreement of these 
data with those of the present analysis demonstrates the internal consistency of the results obtained by 
different methods. 

The sin(/) modulation amplitudes of the azimuthal asymmetry are shown in Fig. [5] Such a modulation is 
expected as a combined effect from the twist-3 PDFs hi and f^ entering doot as well as from the twist-2 
transversity PDF hi and Sivers PDF f^j- entering daor (see Eqs. ^), all contributing to the azimuthal 
asymmetries with a factor Mx/Q. The individual PDF contributions can not be separated within a single 
experiment. The observed x dependence of this amplitude is less pronounced in the COMPASS data 
than in the HERMES |[24l data.^ The latter are obtained for leading pions, while our data include all 
hadrons and cover a much wider range in x, Q^ and W^. When restricting our kinematic region to that of 
HERMES for the amplitude a**'"*^ compatible results are obtained. 

The amplitudes of the sin(20) modulation shown in Fig. [6]ai-e small, consistent with zero within the 
en^ors. They are due to the PDF h-^^ in do^i (see Eqs. ^) which is approximately linked [29] to the 
transversity PDF hi by a relation of the Wandzura-Wilczek type. 

The data on the modulation amplitude a'^'"^'^ shown in Fig. |7]are compatible with zero as our results 



■'The sign of HERMES data was inverted in order to match our definition of spin asymmetry by Eq. (|4j. 
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Fig. 4: Dependence of the 0'=°"'^' parameter on the kinematic variables. The values of A!'. ,,(x) from Ref. Q are 



also shown. 
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Fig.5: Dependence ofthe modulation amplitude a^'"*^ on thekinematic variables and similar data of HERMES 124] 
for identified leading pions. 
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Fig. 6: Dependence of the modulation amplitude a^^'^^'l' on the kinematic variables. The x dependence of this 
amplitude is compared to the data and calculations by H. Avakian et al. l29l for HERMES kinematics^' for positive 
(solid line) and negative (dashed line) hadrons. 
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Fig, 7: Dependence of the modulation amplitude a^^^^^ on the kinematic variables. 
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Fig. 8: Dependence of the modulation amplitude a'^°^'l' on the kinematic variables. 
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Fig. 9: Comparison of results for the h azimuthal asymmetry parameter a'^°"^'(jc)/Do, fl'^°"^'(z) and fl'^°"^'(/9j) for 
different regions: z = 0.05-0.9 and z = 0.2-0.9. 



on the amplitude of the sin(30 — (j)s) modulation extracted from transversely polarised deuterons 11261 
(0s = 0, n for longitudinal target polarisation). This modulation would be due to the pretzelosity PDF 
hjf in dooT and is suppressed by the factor tan dy ~ xM/Q. 

The cos modulation of the azimuthal asymmetries for a longitudinally polarised target is studied here 
for the first time. The data presented in Fig. [8] are consistent with no variations of the modulation 
amplitudes vs. x, z and pj. This amplitude is proportional to the muon beam polarisation and would be 
due to a pure twist-3 PDF gj^ in Aon, an analogue to the Cahn effect QUI in unpolarised SIDIS Il28l . and 
g\T in dOiT, suppressed by tan By ~ xM/Q. 



6 Stability of the results and systematics 

In the above analysis the z cut (z > 0.2) has been applied to reject hadrons originating from the target 
fragmentation region. To check if a lower cut can affect the results presented in Table 1 and in Figs. |4]-[8j 
we have repeated the analysis with the cut z > 0.05 and obtained consistent results. The results for the 
two different z cuts are compared in Fig. |9]for the asymmetry parameters a^°^^^{x)/DQ, a'^°"'^'(z), and 
a!^™^^^{pj). The x dependence of this parameter does not depend on the applied z cut, i.e. there is no 
influence of the target fragmentation region on the x dependence for the selected sample of events down 
to z = 0.05. The observed pj dependence looks smoother with higher statistics. 

The compatibiUty of the final results on the azimuthal asymmetries obtained with the data taken in 2002- 
2004 has been checked by building distributions of "pulls" 



{a) 



(10) 



'{a) 



where a,- is the asymmetry in a given year for a given hadron charge and bin of the kinematic variables x, 
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Fig. 10: Distribution of pulls for measurements of a,-. 



z, and pj. The value {a) is the corresponding weighted mean over three years and o is the statistical error. 
The overall distribution of pulls for all a, measurements (540 entries, i.e. 5 asymmetries for positive and 
negative hadrons, 3 variables with 6 bins and 3 yeai^s of data taking) is shown in Fig. [TOl As expected, 
the distribution of pulls follows a standard distribution with mean close to zero and sigma equal to unity 
indicating that the fluctuations in the data are only statistical. 

The double ratio Eq. ^ has been used to extract the asymmetries because of the cancellation of accep- 
tance and flux. To check this cancellation, two different double ratios have been constructed using the 
same number of events, namely 



and 



^+(0) 



F_(0) 









(11) 



(12) 



for the positive and negative solenoid field orientations, respectively. Substituting into Eqs. (fTTl [T2l) 
the number of events by the corresponding expressions given in Eq. O, one could see that these ratios 
depend neither on the acceptance nor on the DIS cross-section. The above equations should reflect the 
relative integrated muon fluxes squai^ed for the measurements with different target polarisation but the 
same solenoid field orientation and they are a good check of acceptance cancellation in Eq. Q. 

Indeed, the ratios Eq. (fTTl [T2l) for the selected events integrated over all kinematic variables are found 
to be independent of <p and the field orientations and are the same for /i+ and h^ (i.e., 0-dependent 
acceptances CUip) are really cancelled). This is illustrated in Fig. [TT] where the weighted sums of F = 
f+(0)©F_((^) for negative and positive hadrons are shown together with results of their fit by a constant. 

The compatibility and stability tests indicated no systematic effects. It is therefore assumed that the 
systematic uncertainties of the measured asymmetries due to variations of the acceptance are smaller 
than the statistical errors. 

The uncertainties due to the target and beam polarisation are estimated to be of the order of 5% each. 
The uncertainty of the dilution factor, which takes into account the target material composition, is of the 
order of 2% |6]. Each of them introduces the corresponding multiplicative uncertainty in the asymmetry 
measurement. When combined in quadratures, these errors give a global systematic multiplicative un- 
certainty of less than 7%. The errors due to the additive radiative correction are negligible in the used 
kinematical region. 
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Fig. 11: Double ratios of event numbers for the weighted sum F = F+{<p) ®F^{(j)) demonstrating the (j) stabiUty 
of the data taking rates for negative (top), positive (middle) and all hadrons (bottom). The parameters pQ represent 
the results of fits by a constant. 

7 Conclusions 

Azimuthal asymmetries a{<p) were studied in the production of positive and negative hadrons by 160 GeV 
muons scattering off the longitudinally polarised deuterons. Integrated over the variables x, z, and pj, 
all 0-modulation amplitudes of a{<p) are. consistent with zero, while the 0-independent parts of the a{<p) 
distributions differ from zero and are equal for positive and negative hadrons within the statistical errors. 
In the study of a (0) as a function of jc, z, and pj the following results are obtained. The -independent 
terms a^°'^^\x) of a(^) ai^e in agreement with our earlier results Q and other published data on A\^, 
calculated by other methods and using different cuts. The amplitudes a'*'"'^ are small and in general com- 
patible with the HERMES data IH. The amplitudes a''''^1'{x,z,p'j-), a''''^'t'{x,z,p'j-) and a'""^{x,z,pf^) 
are consistent with zero within statistical en^ors. 

These data will be useful to constrain models of the nucleon structure. The present parton model descrip- 
tion of the SIDIS cross-sections involves a considerable number of PDFs depending on the longitudinal 
and transverse components of the nucleon spin. We believe that our data will help to assess which PDFs 
are important in the description of the nucleon structure. 
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